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Abstract. 

Systematic effects on HI and Ha long-slit observations make a mea- 
surement of the inner slope of the dark matter density distribution diffi- 
cult to determine. Halos with constant density cores and ones with r~^ 
profiles both appear consistent with the data, although constant den- 
sity cores generally provide better fits. High-resolution, two-dimensional 
velocity fields remove most of the systematic effects, yet as a result of 
noncircular and random motions the inner slopes still cannot be accu- 
rately measured. Halo concentration parameters provide a more useful 
test of cosmological models because they are more tightly constrained 
by observations. The concentration parameters for LSB galaxies appear 
consistent with, but on the low end of the distribution predicted by CDM. 



1. Introduction 

Studies of dark matter in spiral galaxies have been plagued by uncertainties in 
the mass-to-light ratio (M/L) of the stellar disk, resulting in a large degeneracy 
in the mass modeling (e.g., van Albada et al. 1985). This degeneracy can 
largely be avoided by studying dwarf and low surface brightness (LSB) galaxies. 
Even though the central parts of the observed rotation curves can in principle 
be explained by scaling up the contribution of the stellar disks (Swaters 1999; 
Swaters, Madore, & Trewhella 2000), the inferred stellar M/Ls are much higher 
than expected from population synthesis modeling (e.g.. Bell &: de Jong 2001). 
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For reasonable M/Ls, these galaxies are dominated by dark matter at all radii, 
making dwarf and LSB galaxies ideal for studying the properties of dark matter. 
The power law slopes of the central dark matter density distributions p(r) oc 
are of particular interest. Cosmological simulations indicate a depends on 
the nature of dark matter. A measurement of a may thus provide constraints 
on the nature of dark matter and theories of galaxy formation. Unfortunately, 
it is difficult to measure the inner slope a accurately, both from simulations 
(e.g.. Power et al. 2003; White, this volume; Navarro, this volume), and from 
observations, as is discussed in detail below (see also de Blok, this volume). 

2. Systematic effects on rotation curves 

The main assumption in using rotation curves to infer the properties of dark 
matter halos is that the rotation velocities reflect the gravitational potential. 
However, observational effects, such as poor spatial resolution or slit offsets, and 
intrinsic properties of the galaxies, such as pressure support or noncircular mo- 
tions, may affect the derived rotation curves and lead to systematic uncertainties 
in the derived dark matter properties. 

2.1. HI observations 

Earlier HI studies found that dwarf and LSB galaxies have slowly rising inner 
rotation curves (see e.g., Swaters 1999 and references therein). Because of their 
shallow slopes, these rotation curves were found to be inconsistent with the 
steep halos predicted by cold dark matter (e.g., Navarro et al. 1996; McGaugh 
Sz de Blok 1998). More recent studies raised the concern that some of these HI 
rotation curves may be affected by beam smearing, and found that the rotation 
curves of dwarf and LSB galaxies rise more steeply when beam smearing is taken 
into account (Swaters 1999; Blais-Ouellette et al. 1999; Swaters et al. 2000). 

However, beam smearing in itself does not constitute an unsurmountable 
problem. It is possible to correct the measured rotation curves for the effects of 
beam smearing (Swaters 1999), or its effects can be incorporated in the model- 
ing procedure (e.g., van den Bosch et al. 2000). Nonetheless, because of their 
relatively low resolution, HI observations may not be best suited to determine 
the inner rotation curve slopes. High angular resolution observations, such as 
Ha long slit or Fabry-Perot spectroscopy, thus seem more suited. 

2.2. Ha long-slit observations 

Despite the high spatial resolution of Ha long-slit observations, it remains dif- 
ficult to measure a accurately. Different studies find conflicting results: some 
find LSB galaxies to be consistent with steep inner slopes (e.g., Swaters et al. 
2003a, hereafter SMvdBB), yet others find them inconsistent with steep slopes 
(e.g., de Blok et al. 2001a,b; de Blok & Bosma 2002; Marchesini et al. 2002). 
This inconsistency remains even when identical datascts are used (SMvdBB). 

This apparent discrepancy may well be the result of systematic effects on 
the long-slit observations, which all lead to an underestimate of the inner slope. 
Several effects may contribute: 1) Seeing tends to wash out the inner gradient of 
the rotation curve, leading to an underestimate of the inner slope. 2) If the slit is 
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opened wider, parts of the galaxy with lower radial velocities will also be visible 
through the slit, diluting the inner gradient and resulting in an underestimate of 
the inner slope. 3) Slit offsets (the combined results of telescope pointing errors, 
errors in the coordinates of the optical centers, and offsets between optical and 
dynamical centers) will cause the steepest gradient in the rotation curve to be 
missed, leading to an underestimate of the inner slope. 4) At inclinations close 
to edge-on, one line of sight samples a large part of the disk, resulting in a 
range of observed velocities. The rotation velocity is represented by the extreme 
velocity along the linc-of-sight. Obviously, the usual method of deriving the 
rotation velocity by Gaussian fits or a intensity weighted mean will result in an 
underestimate of the rotation velocity, which in turn leads to an underestimate 
in the inner slope. On top of that, extinction and the distribution of Ha (e.g., a 
ring) may also lead to an underestimate of the inner slope. 5) If the Ha emission 
is concentrated toward the edge of the slit, the derived velocities will be lower, 
especially near the center of the galaxy, leading to an underestimate of the inner 
slope. In summary, all systematic effects lead to an underestimate of the inner 
slope, so there clearly is an strong observational bias against finding steep slopes 
from long-slit observations. 

3. The inner slope: poor constraints on cosmological models 

In the literature, three principal ways of measuring the inner slope of the halo 
density profile have been used: inversion of the rotation curve (e.g., de Blok 
et al. 2001, SMvdBB), fitting a power law to the observed rotation curve (e.g., 
Simon et al. 2003), and fitting a general mass model to the observed rotation 
curve (e.g., van den Bosch et al. 2000, SMvdBB). 

3.1. Rotation curve inversion 

If one assumes the disks are dynamically insignificant, and if one assumes fur- 
thermore that the dark matter has a spherically symmetric distribution, it is 
possible to recover the density distribution of the dark matter from the observed 
rotation curve in a non-parametric way. Prom = A-irGp and $ = —GM/r 
the density distribution is given by 



The inner slope a can be measured by inverting the measured rotation curve 
using Eq. 1, and fitting a simple power-law p oc r"""" to the density distribution 
inside the break radius (the radius where the density profile turns from its inner 
slope to its outer slope). 

In Fig. 1 we present the distribution of measured inner slopes for the galaxies 
in the samples of dc Blok et al. (2001a), dc Blok & Bosma (2002), and SMvdBB). 
Overplotted on the histograms are the expected distributions of measured inner 
slopes based on detailed modeling of the observations, including the effects of 
distance, seeing, slit width, and slit alignment errors. See SMvdBB for more 
details on the modeling. 
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Figure 1. The distribution of measured inner slopes for the galaxies 
in the samples of de Blok et al. (2001a), de Blok & Bosma (2002), 
and SMvdBB), divided into non-barred, barred, and edge-on galaxies, 
both for galaxies within 2000 km s~^ (left), and for galaxies beyond 
2000 km s~^. Dotted lines represent the expected distributions of inner 
slopes for a = 0, the solid lines from left to right give the expected 
distributions for a = 1, with a slit alignment error of 2", 1" and 0". 



For the non-barred galaxies the measured inner slopes span a wide range 
from a = to a = 1. The distribution of a^n for the non-barred galaxies in 
the combined sample appears inconsistent with the distribution expected for 
halos with a = 0. However, it is possible that the stellar disks do contribute 
significantly and this may explain the wing to high am even if the halos have 
a = 0. The models with a = 1 appear to underpredict the number of galaxies 
with am near zero, unless the error in the slit positioning is about 2". An error 
of 2" in the slit positioning may not be unreasonable given the large angular size 
of these galaxies and their irregular, diffuse and low surface brightness nature. 

Although based on only a few galaxies, the barred galaxies on average ap- 
pear to have a somewhat lower am- This may be an indication that the non- 
circular motions associated with the bar have affected the inner slopes. Another 
striking feature in Fig. 1 is the large difference in the distribution of measured 
inner slopes between edge-on galaxies and non-barred galaxies. This difference 
is most likely artificial: the method with which the rotation curves were derived 
lead to a severe underestimate of the inner slopes (SMvdBB). 
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Figure 2. Prom top to bottom the best fit c, best fit V200, and reduced 
X^for the generic halo fits. The the solid line represents X = (for the 
i?-band), the dotted line is X = 0.5, the short dashed lines is X= = 1 
and the long dashed line is X = 2. Halos with < a < 1 are consistent 
with most galaxies, although halos with a = generally provide better 
fits. 



Although the inner slope derived from rotation curve inversion has as ad- 
vantage that it is a parameter-free measurement and that it directly measures 
the inner slope of the halo, it is particularly sensitive to systematic effects men- 
tioned above because the inner slope is usually measured from only the innermost 
points. As a result, it is not possible to measure the inner slopes accurately, and 
our conclusion is that inner slopes in the range < a < 1 are consistent with 
the observations discussed here. 

3.2. Fitting a power law to the rotation curve 

Some authors have determined the inner slope by fitting a power law slope to 
the entire observed rotation curve (e.g., Simon et al. 2003). Although this gives 
a more accurate measurement than rotation curve inversion, it is seems unlikely 
that an inner slope derived from a fit to the entire rotation curve can be directly 
compared to the numerical simulations, because in the simulations the slopes 
changes with radius. A measurement of the true inner slope could be derived 
by fitting a power law slope to the inner parts of the rotation curve only, but 
that would make the measurement sensitive to the systematic effects described 
in the previous subsection. 

3.3. Fitting a general density profile to the rotation curve 

The inner slope can also be measured by fitting a halo model to the rotation 
curve based on a general density profile of the form: 

^^"■^ " (r/r,)"(l + r/r,)3-"' 

where po is the central dark matter density, and a scale radius. This density 
distribution changes from p oc r~" for r <C to p oc for r 3> rg. For a = 
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Figure 3. The shaded histogram gives the distribution of best fit c 
parameters for galaxies in the samples of de Blok et al. (2001a), de 
Blok & Bosnia (2002), and SMvdBB. The arrows indicate the best fit 
c parameters derived from new SparsePak data (Swaters et al. 2003, 
Swaters et al. in preparation). 



the density profile thus has a constant density core and becomes comparable to 
that of the pseudo-isothermal sphere, while for a = 1 it reduces to the NFW 
profile (Navarro et al. 1997). 

The density and scale radius derived from such mass models are less sensi- 
tive to the effects of seeing, slit width and slit offsets, because the entire rotation 
curve is used to determine the parameters, not only the innermost part. Conse- 
quently, these mass models are not very sensitive to the inner slope of the dark 
matter halo, and for most galaxies we find that mass models based on a = 1 
provide fits of similar quality as those based on a = (see examples in Fig. 2). 
This is not only true if the contribution of the stellar component is ignored, but 
also for models in which the stellar component play a modest role. 

From these general mass models we find that for most galaxies the inner 
slope is unconstrained, and fits of similar quality are obtained for the range 
< a < 1, although a = is somewhat preferred. Inner slopes as steep as 
a = 1.5 are clearly ruled out. 



4. Halo concentration: useful constraints on cosmological models 

A different way of testing the predictions of CDM simulations is through the 
halo concentration parameter c. For an NFW halo (Eq. 2 with a = 1), c = 
f20o/fs, where r2oo is the virial radius. An important advantage of using the halo 
concentration parameter is that it is less sensitive to systematic effects because 
it is derived from a fit to the rotation curve as a whole. The expected values 
for c depend on the assumed cosmological parameters. For a ACDM cosmology, 
expected values have a 2a range from 5 to 25, with an average of around 10 
to 15 (Navarro et al. 1997, Bullock et al. 2001). In Fig. 3 the distribution of 
best fit c values is shown; the uncertainty on c from rotation curve fits is about 
20% to 40% for most galaxies. The peak near c = 1 is likely to be artificial, 
resulting predominantly from galaxies in which the fit is unconstrained, such as 
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Figure 4. Radial velocity relative to the systemic velocity, normal- 
ized by rotation velocity, as function of angle in the plane, for DDO 39 
(obtained with SparsePak on the WIYN telescope, Swaters et al. 2003). 
The left panel shows points within r < 1 kpc, the right one points be- 
tween 1.5 < r < 2.5 kpc. The solid line represents circular rotation. 



galaxies with rotation curves that extend less than two disk scale lengths, and 
from edge-on galaxies with possibly incorrect rotation curves. The rest of the 
distribution shows a large spread, and seems somewhat skewed toward smaller 
values of c than expected in ACDM. 

Low values for c do not necessarily indicate an inconsistency with ACDM. 
For example, a bias toward low c values in LSB galaxies could be explained if 
LSB galaxies preferentially form in low density halos. In addition, the concen- 
tration parameter also depends on the slope of the power spectrum of density 
fluctuations, and c ~ 5 is in agreement with models in which structure formation 
on small scales is suppressed (Zentner & Bullock 2002). 

5. High Resolution, Two-Dimensional Velocity Fields 

To avoid systematic eff'ects of long-slit and HI observations, and to map possible 
non-circular motions, high spatial resolution, two-dimensional velocity fields are 
needed. Several such studies are now available (e.g., Blais-Ouellette et al. 2001; 
Bolatto et al. 2002; Swaters et al. 2003b; Simon et al. 2003). 

Swaters et al. (2003b) studied the kinematics in the core of the LSB galaxy 
DDO 39 with the SparsePak integral field spectrograph on the WIYN telescope 
(Bershady et al., submitted), and found random motions with an average ampli- 
tude of about 8 km s~^ throughout the galaxy, and in addition detected distinct 
noncircular motions as can be seen in Fig. 4. We (Swaters et al. in preparation) 
have obtained data for an additional 6 LSB galaxies with SparsePak. Prelim- 
inary analysis of these data again reveal significant random and noncircular 
motions. The origin of either component is unclear, but they cause inner slopes 
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in the range < a < 1 to be consistent with the data. Therefore, even though 
these data are a significant improvement on previous HI or Ha long-slit data, the 
inner rotation curve, and hence the inner slope of the dark matter distribution, 
is still uncertain due to significant noncircular and random motions. 

Fortunately, the availability of two-dimensional velocity fields makes it pos- 
sible to determine the kinematic center and estimate the influence of noncircular 
motions on the derived rotation curves, allowing for an accurate measurement 
of all but the inner rotation curve. Consequently, halo parameters, such as c, 
that depend on the entire rotation curve shape can be determined accurately 
(with a typical uncertainty on c of about 10% to 20%). For DDO 39 and for the 
new SparscPak data, we found c values that are on the low side when compared 
to simulations (see Fig. 4). 

6. Conclusions 

The inner slope a is difficult to measure from HI and Ha long-slit observations 
because of systematic effects. High-resolution, two-dimensional velocity fields 
provide a major step forward, but deviations from circular motions still make 
an accurate measurement of a difficult. A range of inner slopes < a < 1 
appear consistent with the data, although models with a = generally pro- 
vide somewhat better fits. Thus, it seems that an accurate measurement of 
a and, from that, a tight limit on the nature of dark matter may be difficult 
to achieve observationally from Ha emission, which is susceptible to kinematic 
perturbations. 

Halo parameters that have been derived from the entire rotation curve in- 
stead of only the central parts, such as the halo concentration parameter c, may 
provide a more robust test for cosmological models. In the case of NFW ha- 
los, the halo concentration parameters of LSB galaxies tend to be lower than 
expected from cosmological simulations. 
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